Disordered, Spin Liquid and Valence-Bond Ordered Phases of the Kagome Lattice 
Quantum Ising Models With Transverse Field and XXZ Dynamics 
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General conditions in which disordered, spin liquid, and valence-bond ordered phases occur in 
quantum Ising antiferromagnets are studied using the prototype Kagome lattice spin models. A 
range of quantum dynamical processes in the Ising model, with and without total Ising spin con- 
served, are analytically shown to yield all three characteristic quantum paramagnetic phases in the 
Kagome system. Special emphasis is given to the XXZ model that can be sensibly compared to 
the Kagome lattice Heisenberg antiferromagnet. It is explicitly demonstrated that the total-spin- 
conserving dynamics can yield a resonant valence bond (RVB) liquid phase with very short-ranged 
correlations, but also a valence-bond ordered phase compatible with the one proposed to explain the 
seemingly gapless singlet states of the Heisenberg antiferromagnet on the Kagome lattice. Likely 
consequences for generic spin models are discussed. The analysis combines compact U(l) gauge 
theory, duality transformations, lattice-field-theoretical methods, and variational approach. 
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I. INTRODUCTION 

The quest for a quantum spin liquid has been a ma- 
jor pursuit in condensed matter physics ever since the 
Anderson's proposal for its existence and relevance to 
the cupratesi Since then, some new materials with geo- 
metric frustration emerged as promising candidates for 
the spin liquid physics at low temperatures. Among 
two-dimensional systems the Mott insulators CS2CUCI4 
(Ref^i) and k - (BEDT-TTF)2Cu2(CN)3 (Ref.^), based 
on the triangular lattice, experimentally exhibit uncon- 
ventional magnetic behavior in certain circumstances, 
without detectable symmetry breaking. Numerical stud- 
ies also reveal possible spin liquid in realistic systems, 
ranging from a Wigner crystal near melting^ to various 
antiferromagnets with multiple-spin exchange^*^ How- 
ever, there is still no unambiguous experimental evidence 
that a spin liquid is found in any of these or similar cases. 
Even theoretically, much more is known about general 
properties of the spin Iiquid£*2i2ii2iii then about funda- 
mental and microscopic circumstances needed for its re- 
alization. Gaining more microscopic insight can only be 
useful for understanding the unconventional magnetic be- 
havior found in various experiments, and perhaps also for 
clarifying available options for quantum computing^Sii^ 

In this paper we attempt to learn general lessons on 
the role played by the lattice structure, symmetries and 
type of dynamics in shaping the phases of the frustrated 
quantum magnets. The present analysis is focused on the 
prototype Kagome lattice Ising antiferromagnets, but the 
range of explored models admits all characteristic quan- 
tum paramagnetic phases: disordered, spin liquid, and 
valence-bond solid. Through connecting these outcomes 
with both the fundamental and microscopic properties 
of the models, and with information that emerges from 
calculations, we can deduce some conclusions of broader 
significance for the spin models on other lattices. 

The Kagome lattice is an excellent choice for this pur- 
suit because it is one of the few simple spin systems (with 
only nearest-neighbor interactions) where disordered and 



spin liquid phases are believed to exist. On the other fre- 
quently studied lattices it usually takes further-neighbor 
and multiple-spin exchange to destabilize the zero tem- 
perature Neel order. The situation is somewhat bet- 
ter when the frustrated Ising antiferromagnets are con- 
cerned, but even then the quantum dynamics typically 
leads to a paramagnetic ground state that breaks lattice 
symmetries (such as on the triangular and fully frustrated 
square lattices). 

The Kagome antiferromagnets are unique among 
two-dimensional spin systems in that they not only 
exhibit a promising spin liquid- like behavior, but 
also sometimes hold additional surprises. The ex- 
perimental research on the spin = I Ka gome 
layered materials SrCr9pGai2-9pOi9 (SCGO )^"^!^^ and 
Ba2Sn2ZnGa3Cr7022 (QS ferrite)i& discovered a spin- 
glass together with a heat capacity that is not thermally 
activated, and largely not dependent on the magnetic 
field at low temperatures. The numerical exact diag- 
onalization studies of the S = j Heisenberg antiferro- 
magnet on small Kagome samples^^'^'* revealed a dis- 
ordered ground state, and a seemingly gapless band of 
numerous singlet excitation that fill the spectrum below 
the finite spin-gap. In comparison to the other quan- 
tum paramagnets^SiSS this one appeared qualitatively 
different, and called for classification as a new and ex- 
otic kind of spin liquid. Such kind of spectrum, in which 
there seems to be no gap in a completely disordered 
phase, is still not understood. The physical pictures 
proposed by various theoretical efforts favored a spin 
liquid but some new ideas open up a pos- 
sibility of a valence bond crystal with very large unit-cell 
that gives rise to an extremely low energy scale for the 
singlet degrees of freedomi^SiSl 

The transverse field Ising model on the Kagome lat- 
tice is another example of how special this lattice is. The 
Monte-Carlo simulations^^'^^ pointed out that unlike the 
other Ising systems, Kagome prefers not to order even 
for small transverse fields. One perspective in explain- 
ing the reason for this has been taken in the Ref.'^*' In 
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FIG. 1: Kagome lattice is a corner-sharing two dimensional 
lattice. The frustrated units are triangular plaquettes, and 
they are only minimally connected into the lattice. 

this paper, a different perspective will be given. The 
present approach is also extended to the XXZ model, 
which can be regarded as the Heisenberg model with 
easy-axis anisotropy. 

II. MODELS AND OVERVIEW 

In this paper we analyze the nearest-neighbor spin S = 
^ quantum Ising antiferromagnets on the Kagome lattice 
(Fig. Two kinds of spin dynamics will be explored, 
represented by the following simple models: 

• transverse field Ising model (TFIM): 

i? = J.^5f^|-r^5f ; (1) 

iv) i 

• Heisenberg model with easy-axis anisotropy (XXZ) : 

In contrast to the transverse field case, the XXZ dy- 
namics preserves total Ising magnetization, making the 
Hamiltonian (j^J symmetric under global spin-flip. Fur- 
thermore, the transverse field gives rise to the most local 
kind of spin dynamics, while the XXZ dynamics involves 
pairs of spins, and thus introduces some correlation. It 
will become apparent that these two fundamental differ- 
ences yield very different low energy physics. The conse- 
quent analysis will also admit introduction of other dy- 
namical processes, spatially extended to larger clusters 
of spins, but consistent with the symmetries of these two 
basic models. 

The calculations in this paper are restricted to 
weak dynamical perturbations of the pure Ising model: 
r, Jj_ <C Jj. This limit is a combination of analytical con- 
venience, and essential physical interest in the context of 



frustrated magnetism. The main question being asked is 
how the quantum fluctuations (created by weak dynam- 
ical perturbations) lift the degeneracy of the pure Ising 
model. Is the ground state ordered like in many other 
Ising systems (order-by-disorder on triangular and fully 
frustrated square lattices)? Under what general circum- 
stances a completely disordered ground state is possible, 
with or without non-trivial topology? In attempt to an- 
swer these questions, we will formulate a lattice field the- 
ory and apply to it a technique specialized for frustrated 
systems, but otherwise analogous to the usual mean-field 
approach in the unfrustrated problems. Namely, instead 
of finding the mean-field solutions that minimize energy 
in various parameter regimes, we will seek solutions that 
maximize "entropy" of quantum fluctuations. -21 When 
needed, those solutions will be subject to a verification 
of stability. This will provide a reliable picture of the 
phases that exist in our models. 

Physics of the TFIM is trivial when the transverse 
field r is large, while in the limit F ^ the quan- 
tum dynamics, as a matter of principle, has a chance 
to yield interesting valence bond ordered or disordered 
ground states after lifting the huge degeneracy of the 
pure Ising model. Even though this issue has been al- 
ready understood for the Kagome TFIM,'^" the following 
approach is going to bring some more insight: it will al- 
low us to propose certain variational wavefunctions. Of 
course, much closer to the true challenge posed by the 
Kagome Heisenberg antiferromagnet is the XXZ model, 
which can be regarded as its anisotropic version. In the 
limit J_L <^ Jz that we are going to explore the easy-axis 
anisotropy is strong, but sensible comparison with the 
isotropic case will become apparent. In particular, it will 
be shown by variational arguments that the short-range 
valence bond states arise most naturally in the easy-axis 
XXZ models on the corner-sharing lattices, which to- 
gether with the lattice-field-theoretical calculations will 
lead to the valence bond ordered state already proposed 
for the isotropic Kagome antiferromagnet^ A spin liquid 
phase will also appear as a possibility in closely related 
models. 

This paper is organized in two main sections, each di- 
vided into several subsections. The TFIM model is dis- 
cussed in the section IIIII while the XXZ model is studied 
in the section IIVI Initial discussion of the XXZ model re- 
lies heavily on the definitions and ideas introduced in the 
TFIM section (sections IIIII through IIIID|I . Readers in- 
terested only in the physical nature of the valence-bond 
ordered and spin liquid phases of the XXZ models can 
skip all calculations and go directly to the section IIV Fl 
All conclusions arc summarized and bigger perspective is 
taken in the Discussion. 
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III. TRANSVERSE FIELD ISING MODEL 

We start from the Kagome lattice Ising model in a 
weak transverse field F ^ J^: 



(3) 



Let us first understand the ground states of the pure Ising 
Hamiltonian (F = 0). They are the least frustrated states 
in which the number of frustrated bonds (two aligned 
spins) is minimized. If every frustrated bond is visual- 
ized by a dimer, then every appropriate dimer covering 
determines a spin configuration up to a global spin flip. 
Consider a loop on the Kagome lattice (Fig. ^ . The un- 
frustrated bonds on the loop mark locations where the 
two neighboring spins on the loop have different orienta- 
tion. When going one full circle around the loop one ends 
at the same spin from which one started, so that the num- 
ber of times the spin orientation is changed must be even. 
Therefore, every loop contains an even number of unfrus- 
trated bonds, and the parity of the number of dimers on 
the loop is determined by the loop size. The number of 
dimers on the triangular (hexagonal) Kagome plaquet- 
tes must be odd (even). This is the only constraint for 
the Kagome lattice dimer coverings that correspond to 
arbitrary spin states. 

It will be convenient to immediately switch to the dual 
picture. Duality between the Kagome and dice lattices is 
depicted in the Fig. |31 Since every Kagome bond corre- 
sponds to one dice lattice bond, the frustrated bonds can 
be represented by dimers on either lattice. An example 
is shown in the Fig. 21 There must be an odd (even) 
number of dimers emanating from every 3-coordinated 
(6-coordinatcd) dice lattice site. The number of dimers 
(and thus frustration) is minimized if there is exactly 
one dimer emanating from every 3-coordinated dice lat- 
tice site. This condition flxes the number of dimers in 
the least frustrated states, since the dice lattice is bi- 
partite. The degeneracy of the least frustrated states is 
apparently huge. 

Small dynamical perturbations will mix the least frus- 
trated states and lift their immense degeneracy. In prin- 
ciple, one can perturbatively derive an effective theory 




FIG. 2: Every loop holds an even number of unfrustrated 
bonds. Dimers represent the frustrated bonds, or pairs of 
aligned spins. 




FIG. 3: Duality between the Kagome (solid line) and dice 
(dashed) lattices. Duality transforms a Kagome site into the 
dice plaquette inside which it sits, and vice versa. The 3- 
coordinated dice sites and Kagome triangles transform into 
each other, as well as the 6-coordinated dice sites and Kagome 
hexagons. Every Kagome bond intersects one dual dice bond. 




FIG. 4: Frustrated bonds represented by dimers on the dice 
lattice. 



that describes dynamics at energy scales well below J^. 
This effective theory lives in the Hilbert space spanned 
only by the least frustrated states. Therefore it takes 
form of a soft-core quantum dimer model on the dice 
lattice, where exactly one dimer emanates from every 3- 
coordinated dice site, while an arbitrary even number of 
dimers emanate from every 6-coordinated site. For our 
purposes it will be sufRcient to concentrate just to the 
first order of degenerate perturbation theory: 



H, 



eff ■ 



2^i\<'){> 



+ h.c. 



h.c. 



(4) 



The dimer dynamics consists of two different flips on the 
dice plaquettes that are consistent with minimum frus- 
tration (see Fig. inj. Note that these two processes in- 
volve flipping of only one spin on the Kagome lattice. 
For the purposes of simplicity and staying close to the 
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P,<1- 

Pi 

P6 

{ij) 



Kagome lattice sites, or dual dice plaquettes 
dual dice lattice sites, or Kagome plaquettes 
3-coordinated dice sites, or Kagome triangles 
6-coordinated dice sites, or Kagome hexagons 
Kagome lattice bonds 
dual dice lattice bonds 



TABLE I: Notation for the Kagome and dual dice lattices 



original spin dynamics, we will not consider a more gen- 
eral dimer model with different energy scales for the two 
types of flips. Since the dice lattice is bipartite, it is 
possible to apply the standard techniques and cast this 
dimer model as a compact U(l) gauge theoryi^ From 
that point on, duality transformations and lattice-field- 
theoretical methods are at disposal to study the possible 
phases. After understanding the phases that the lattice 
theory yields, it will be apparent that sufficiently small 
higher order perturbations cannot destabilize them. 

Some comments are in order before proceeding. First, 
the dimer representation is insensitive to the global spin 
flip. This is of no concern for the models studied here, 
since magnetically ordered phases will not be found. Sec- 
ondly, the dice lattice dimer model in its own right has 
non-trivial topology, just like the other quantum dimer 
models. However, only one of its topological sectors cor- 
responds to the physical spin states on a torus. There- 
fore, the original spin model need not have any interesting 
topology. 




FIG. 5; The reference bond orientations of the Kagome and 
dice lattices. The bipartite dice lattice bonds are oriented 
from the 6-coordinated site to the 3-coordinated site. Ev- 
ery Kagome bond orientation is locked to the orientation of 
the dual dice bond by the "right hand rule". Note that the 
Kagome orientations circulate around triangles and hexagons 
in the different directions. 



cable to both the dice and Kagome lattices. The bond 
orientation will transform by duality according to the 
"right hand rule". Since the dice lattice is bipartite, we 
will orient its bonds in a natural way, and chose the orien- 
tation to be from the 6-coordinated to the 3-coordinated 
site on every bond. This fixes orientation of the Kagome 
bonds as well, and we show both in the Fig. |S1 

Now we can define the electric field Epq as a vector 
corresponding to the scalar bond energy i?(pq) : 



A. Compact U(l) Gauge Theory 

Calculations in this paper rely heavily on duality be- 
tween the Kagome and dice lattices (see Fig.|2)l. In order 
to facilitate mathematical manipulations, we will treat 
both lattices on the same footing, and regard the pairs 
of objects related by duality as identical. The notation 
that we will use from now on is summarized in the Table 
n Note that, according to this principle, any quantity 
that lives on a Kagome bond equivalently lives on the 
dual dice bond, and may be labeled by either Kagome, 
or dice bond labels. Also, we will apply the following 
convention: if an equation shows a relationship between 
expressions defined on different lattices, the dual lattice 
objects are always implied (for example, Kagome site 
dual dice plaquette, Kagome bond o dual dice bond. . . ). 

Let us also introduce a vector notation. We will dis- 
tinguish vectors Rpq from the corresponding bond scalars 
R{pq) in that the vectors will change sign if the bond ori- 
entation is reversed: Rpq = —Rqp, while the scalars will 
not: R(pq) — R{qp)- In order to establish a formal con- 
nection between the bond vectors and scalars, we assign 
orientation to the lattice bonds. Let the vector rjpq equal 
if the bond (pq) is oriented from p to and —1 other- 
wise. Then, the vectors and corresponding bond scalars 
are related by Rpq = Vpq-I^ipq) ■ These relations are appli- 



^ _ r , vacancy (unfrustratcd bond (pq)) 
^^''^ \ 1 , dimer (frustrated bond {pq)) 

Epq = llpqEl^pq) . (5) 

The Hilbert space of the least frustrated states has re- 
strictions that are easily expressed in the form of the 
Gauss' Law. The number of dimers Epq — rjpq ema- 
nating from any 3-coordinated site is one, and from any 
6-coordinated site an even number (2np^). We use the 
convention that every dice bond is oriented from the 6- 
coordinated to the 3-coordinated site, and write: 

(Vpe) ^Ep^q = '2np^ . (6) 

The interpretation of this Gauss' Law is that there is a 
fixed background charge —1 on every 3-coordinated dice 
site, and a number ^ Up^ ^ 3 of charge 2 bosons on 
every 6-coordinated site. The charged bosons are inde- 
pendent degrees of freedom living on the 6-coordinated 
dice sites. Formally, they emerge because the dice lattice 
dimer model is not hard-core. 
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(a) 



^3 



S3 



q6 < — > P6 



(b) 



■^3 



S3 



FIG. 6: Elementary processes on a dice plaquette that pre- 
serve the minimum frustration, (a) The number of dimers 
emanating from every site is preserved; (b) a pair of dimers 
is exchanged between two 6-coordinated sites. 



Dynamics of the fields and particles can be easily for- 
mulated if the Hilbert space is expanded to allow arbi- 
trary integer strength of the electric field, and arbitrary 
particle occupation. Promoting Epq and into free in- 
tegers between —00 and -f 00 makes it easy to write the 
creation and annihilation operators: exp(±j^pq) for the 
field lines, and exp(±i</3pg) for the particles. The vector 
potential Apq and the boson phase ipp^ are conjugate an- 
gle operators to the electric field Epq and particle number 
ripg respectively: 



[■Apqj-l^pq] — [f 



P6 I "PeJ 



(7) 



After the Hilbert space has been expanded, we must at 
least introduce a large energy cost to all "unphysical" 
states, so that the low energy physics will still correspond 
to the dimer model This is achieved in the large U 
limit through the following new term in the Hamiltonian: 



Hu = Uj2(E^pq^-^y = uJ2e^^ 



const. 



(8) 



(P9> 



(P9> 



The term linear in electric field is a global constant, since 
it expresses the fixed total number of dimers on the dice 
lattice (in the least frustrated states). 

Now we formulate dynamics of Q in the U(l) lan- 
guage. The two processes of interest arc shown in the 
Fig. El Recall that every dimer means i?(pq) = 1, and 
every vacancy — 0, and that ripg is the number of 

dimer pairs emanating from a 6-coordinated site. There- 
fore, we can easily exploit the dice bond orientations, and 
arrange the creation and annihilation operators to de- 
scribe the allowed dimer flip processes (a) and (b) shown 



in the Fig. |B1 

(a) - exp(i^(pg^3))exp(-i^(qg^3)) x 

X exp(j^(,gr3))exp(-iy^(pg^3)) + h.c. 



(9) 



= 2cos(^ 



'P6S3 I 



A 



An 



+ Ar^pg ) 



(6) - exp{iipqj exp{-iipp^) x 

xexp(-i^(p,s3))exp(i^(,,s3>) x 
X exp(i^(,gr3>) exp(-z^(pgr3>) + h.c. 



(10) 



2 C0s((pqg (^pg Apgss -^S3i36 ~^ ■^qers ~^ ■^rspe) 



2 cos I ipq^ - ipp^ + r]p^q^ £{pq)^pq 
pq 



In the last lines of these expressions the sums are taken 
around a plaquette in the counter-clockwise sense; this 
is the lattice circulation, or the curl. The expression 
(|^ is the usual "magnetic" energy, while the expression 
(|10|l is the boson hopping between the neighboring 6- 
coordinated sites. In this paper we will not use a more 
conventional form of particle hopping that involves an 
"integral" of the vector potential along only one path 
between the two sites. In the equation H1U|) we have in- 
troduced two new symbols: e^p,) and r/p^q^. The former 
is needed to correct the signs of Apq that appear in the 
circulation. Note that the signs have been altered with 
respect to the ordinary circulation in (PJ only on one of 
the two paths that connect the two 6-coordinated sites. 
This allows us to choose e^p^^ as shown in the Fig. [7| 
The other symbol, rip^q^ is needed to ensure that the ex- 
pression inside the cosine of H1U|) transforms properly 
when the sites pe and qe are exchanged. This is a new 
vector, defined on the triangular lattice formed by the 6- 
coordinated dice sites, or equivalently the centers of the 
Kagome hexagons. Since it takes the values ±1, it defines 
the bond orientations shown in the Fig. |S| Notice that 
77pggg must be related to £(pq)'. if one takes a closer look 
at the cosines in (IIUII . one can sec that when the boson 
hops from pq to gg 1 the circulation starting from pg must 
first go through e^p^^ = — 1 bonds. When this is satisfied 
for the counter-clockwise circulation, ?7pgqg should be -1-1, 
otherwise it should be —1. Exactly this is achieved by 
relating 77 and e vectors as shown in the Fig. |S1 

Finally, we can summarize the compact U(l) gauge 
theory on the dice lattice. The Hamiltonian is: 



H = uY^E'pq-rY^ 



(pi) 




+ cos ( ipq^ - ipp^ + r]p^q^ '^{pq)Aq 1 

\ PI I - 



, (11) 



FIG. 7: Special bond signs on the Kagome and dice lattices. 
ei^ij) = £(p5) take the value —1 on the emphasized bonds, and 
+ 1 on all other bonds. 




and the Hilbert space is constrained by the Gauss' Law: 



(a) 



(Vp3) E ^P3, = -1 



qep3 



(VP6) E ^' 



Peg — 2r7,pg + 2 



(12) 



For convenience that will become apparent later, we have 
shifted ripg by one in the bottom expression (this sets to 
zero the total background charge on the lattice). In the 
limit of J7 — > oo this Hamiltonian is an exact rewrit- 
ing of the effective theory For finite and large U, 
one perturbatively obtains a theory as a T/U expansion 
that introduces dimer flip processes on larger loops. This 
physically corresponds to further-neighbor and multiple- 
spin exchange that would be also generated by not so 
small r in the original spin model Owing to this 
approximate correspondence between finite U and larger 
r, we might have means to qualitatively see some trends 
beyond very small T/Jz. 



B. Lattice Field Theory 

The path-integral corresponding to Hll|l describes 
a (2-fl)D electrodynamics. All fluctuations are con- 
strained by H12|l . The action will contain a usual Berry's 
phase (we will omit the time index): 



T \(pq) Pe / 



and a potential energy part: 



5pot = C/^rXX^i 



(13) 



(14) 



(P9> 



where 6t is the imaginary time increment. The kinetic 
energy, which involves the cosines in Hll|) , can be brought 




(b) 



FIG. 8: Bond orientations r?pgqg of the triangular lattice 
formed by the 6-coordinated dice sites, or equivalently the 
Kagome hexagon centers. The dice and Kagome bonds with 



'(pg) 



-1 are emphasized. Notice in (b) that every 



triangular lattice bond contains one Kagome site. 



to a more tractable form by applying the Villain's ap- 
proximation. Two new fields will appear and play a sig- 
nificant role: the magnetic field scalar Bi that lives on 
the Kagome sites dual to the dice plaquettes, and the 
particle current jp^qg that lives as a vector on the tri- 
angular lattice bonds. Both will be integer valued, re- 
flecting the compactness of the U(l) gauge theory, and 
the fluctuations of both will be suppressed by the scale 
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I log(r5r/2)|. They take part in the action as follows: 



(15) 




(peqe) 



\ pq 



{pq)-^pq 



After writing this, the angles Apq and ippg can be formally 
integrated out. Fluctuations of the boson phases ipp^ will 
give rise to the particle current conservation law: 



(16) 



Fluctuations of the vector potential will give rise to the 
Maxwell's equation for the magnetic field curl: Apq is 
coupled to the magnetic field Bi, and current jpgge 
(|15|l . as well as the time derivative of the electric field 
Epq in (|13|l . which is the "displacement" current. How- 
ever, this equation will take an unusual form, because the 
particle and displacement currents formally live on dif- 
ferent lattices. The easyest way to derive it is to rewrite 
the terms in which Apq appears using the Kagome lattice 
notation. For this purpose, let us note that the particle 
current jp^qg is related to the triangular lattice bond vari- 
able j(pgqg) by: jp^qg = r]pgqj{peqg), which in turn can be 
regarded as actually living on the sites of the Kagome 
lattice (see Fig.|H}. Therefore, we can label j{pgqg) as ji, 
where i is the Kagome site that sits on the triangular 
bond {peQe)- From ((T^ and ((TB|l we have: 



ApQ^T-Epn — Aij^T-^i' 



Bi ^ ^ Apq — Bi ^ ^ Ai^ 



(17) 



pq 



jpeqeVpeqe ^{pq)'^q — ji ^{ii>-^i 

pq j&i 



The vector potential fluctuations set to zero the sum of 
everything coupled to Aij on every Kagome bond (ij). 
This is the Maxwell's equation: 



(V(ij)) ArE.j = B, - Bj + (j, - J,) 



(18) 



Once the phase and vector potential fluctuations have 
been integrated out, the remaining action contains only 
integer valued fields: 



(19) 



whose fluctuations are subject to the constraints (|12|l . 
(|15|l . and H18|) . For convenience, we have chosen the 



imaginary time increment 5t such that g = USt; this 
will not affect the low energy physics, at least for large 
U. 

Now we can proceed by solving the constraints. To this 
end, we want to completely switch back to the Kagome 
lattice. The duality transformation that follows had been 
worked out in two typical cases, (a) If the compact U(l) 
gauge theory contains no charged particles, then the dual 
theory is an integer-valued height model, (b) If the par- 
ticles and the electric field live on the same lattice, then 
the dual theory has a non-compact U(l) gauge structure, 
and a charged matter field.^i The case (a) emerges from 
the hard-core dimer models on bipartite lattices, while 
the case (b) has been proposed as an approximate de- 
scription of the hard-core dimer models on non-bipartite 
lattices. Our case is somewhat in between. It turns out 
that the dual theory for our case resembles the height 
model. 

It is convenient to redefine the magnetic field and cur- 
rent as time derivatives of two integer-valued "height" 
fields, Xi and Xii 



B, = ArXt 



ji — At-A^ 



(20) 



Now, by introducing the Xp^qg vector analogous to jp^qg, 
we can write solutions of the current conservation Hlt)|) 
and Maxwell's equation lfTF|l : 



'P6 



Xt ~ Xj + (Ai -Xj)+ 



(21) 



(22) 



where Qj is an integer that does not vary with time, and 
will be determined by substituting this expression into 
the Gauss' Law (I12|l . For consistency, let us first rewrite 
the Gauss' Law using the Kagome lattice labels: 



(VOpJ ^ = -2np, - 2 . 



(23) 



The electric field divergence on the 3 and 6-coordinated 
dice sites in (|12ll transforms by duality into the lattice 
curl on the Kagome triangles and hexagons respectively. 
Taking curls of (|22|) will annihilate Xi on all Kagome 
plaquettes, as well as Xi on the Kagome triangles, since 
e^y) is fixed on every triangle (see Fig. [Tj). However, 
Xi will not be annihilated by the curls on the Kagome 
hexagons. One can easily show that: 



E 



'{ij) 



(A. -A,) 



' E/ ^Peqe 
qe&pe 



(24) 
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FIG. 9: One characteristic and periodic configuration for ^(^j) . 
Every dimer represents ^^ij^ = 1, and vacancy (^^ij^ = 0. 



which in turn is equal to — 2npg according to (|21|) . Con- 
sequently, the equations that must satisfy are: 



(VAp3) -1, 



(25) 



There are many possible choices. In order to reveal them 
more tractably, let us use the Kagome bond orientations 
in Fig. 13 and switch to the appropriate bond scalars 
(notice that the counter-clockwise circulations coincide 
with the bond orientations on the triangles, but not on 
the hexagons): 



(26) 



If we decide to use only the values and 1 for , and 
visualize the value 1 as a dimer, we see that every trian- 
gular plaquette must hold one dimer, and every hexag- 
onal plaquette two dimers. One such configuration is 
depicted in the Fig. |5| The other possible configurations 
need not be periodic on the lattice, but "breaking" of the 
translational symmetry is unavoidable as long as are 
integers. 

This concludes the solution of all the constraints (|12|l . 
and (|18|l . The final lattice field theory describes 
fiuctuations of the two integer-valued height fields on the 
(2+l)D Kagome lattice. We obtain the action by substi- 



tuting idni) and ^ into ((B|l : 

S ^ 9j2[E{(^^^^)' + (^r\^)') (27) 

T i 

iv) 

C. Important Properties 

Before proceeding with analysis of fluctuations in the 
lattice theory (|77|l . we have to reveal several of its impor- 
tant properties. We begin by finding all configurations of 



integer-valued Xi = xf'^ and Xi = X\^'' that minimize the 
action. Clearly, xf'^ and A^'^'' should have no time depen- 
dence. However, in terms of spatial variations, there will 
be a large degeneracy. Let us define: 



i(0) 



(0) 



xf +£fe)(Ar-A)"0 + C 



(0) 



(28) 



At a saddle-point, the action reduces to the sum of S,fj 
on all Kagome bonds. The constraints that ^ij must 
obey can be extracted by tracing back to Eij. Re- 
call that this lattice field theory describes a particular 
soft-core dimer model in which — 1 represents a 

dimer. Therefore, the action will be minimized for all 
allowed dimer coverings of the Kagome lattice, by taking 
^(u) ~ ^ ^'^^ ^ dimer, and for a vacancy. There will 
be one dimer on every Kagome triangle, and an arbi- 
trary even number of dimers on every Kagome hexagon. 
Note that every action minimum corresponds to two least 
frustrated spin configurations of the Kagome Ising model 
©, which are related to each other by the global spin flip. 

Shifting the height flelds by x-°^ and A^"^ allows us to 
study fluctuations about a particular saddle-point. The 
action takes a more general form; 



(29) 



r i 



After "summation by parts" , we can write it in a matrix 
form (up to a constant): 



- = E E [2X. E + E ^i^^)^^^ (30) 
+Xi(6Xt - {Xi,r+i + Aj,r-i) -E^-^J' +^(y>^j) 

= x'Cx + ix'i + ex) , 
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where the vectors x ^-nd ^ are arranged as: 



(31) 



I = 



It is crucially important to understand the properties 
of the saddle-point vectors ^ and the coupling matrix C. 
They follow in a straight-forward manner from (|30|l and 
H31|l . but due to tediousncss of algebra, we dcffcr deriva- 
tion to the Appendix Here we will only summarize 
the results: 



• all saddle-point vectors ^ have the same norm: 

const. ; (32) 



• all saddle-point vectors ^ are degenerate eigenvec- 
tors of the coupling matrix C : 



(33) 



In fact, the coupling matrix C is completely dispersion- 
less: its eigenvalues have only frequency dependence, and 
no dependence on spatial wavevectors. If the height fields 
were real instead of integer numbers, there would be six 
localized bare modes per Kagome lattice unit-cell. Two 
would be "gapless" (zero eigenvalue at zero frequency) 
and degenerate, while the other four would be "gapped" 
and degenerate. The "gapless" bare modes are actually 
unphysical, and merely a redundancy of the state rep- 
resentation in terms of the height fields. Exciting them 
at arbitrary places and frequencies does not affect at all 
the only physical quantity in the problem, the bond en- 
ergy Eij (see Fig. I1U|I . Their existence can also be con- 
firmed by counting arguments: there are six bond vari- 
ables (Eij), and two constraints on the Kagome triangles 
(Gauss' Law) per unit-cell, leaving only four independent 
variables per unit-cell. Also note that because of H33|) 
the gapless bare modes do not couple to the saddle-point 
vectors ^ in In other words, they behave much like 

some "gauge" degrees of freedom. 



D. Effect of Fluctuations 

The lattice theory (|29|) that we want to analyze is 
very similar to the simple and well understood integer- 
valued height model. In 2 -f 1 dimensions the height 
model is known to order and give a "smooth" phase. 
This means that fluctuations of the height field Xi are 
such that the average {{Xi~Xj)^) does not diverge as the 
sites i and j go far apart. The question that we now ask 
is whether the Kagome double-height theory (|29|l also 
lives in a "smooth" phase, and what kind of lattice sym- 
metry breaking (if any) is obtained. We will try to find 
answers in a fashion inspired by the Reference i'^* First 




(a) 



(b) 



FIG. 10: Unphysical fluctuations in the lattice field theory 
12711 . The pairs of numbers in the figure indicate all non-zero 
values of (Xii-^O that constitute a redundant integer- valued 
field configuration (up to a multiplicative integer constant). 
One can easily see from 1221 that arbitrary superposition of 
these configurations does not affect the bond energy Eij on 
any bond, ei^ij) is —1 on the bonds of the shaded triangles, 
and -f 1 on the other bonds. 



we investigate which microstates are most frequently vis- 
ited by the system as it fluctuates. This is not a trivial 
problem for the frustrated magnets, since the microstates 
that minimize the action are macroscopically degenerate 
and only the entropy of quantum fluctuations can dis- 
tinguish them. If we somehow found that the preferred 
microstates possessed long-range order, we would have 
to verify whether that order is truly stable against fluc- 
tuations. Note that the formal lack of dispersion in the 
coupling matrix C does not automatically signal localiza- 
tion, and hence the "rough" phase: the bare modes are 
strongly interacting in order to give only integer-valued 
height fields on all sites. 

Finding the microstates that are most frequently vis- 
ited by the system is an extension of the standard mean- 
field approach, where such states simply minimize the 
energy. The appropriate quantity to minimize in general 
circumstances is a "free energy". Physically, when the 
ground state is a superposition of many different states 
(microstates) , the ones with large probability amplitudes 
squared will be characterized by small "free energy" . The 
probability amplitude of a state ip can be in principle cal- 
culated from the imaginary-time (or finite temperature) 
path-integral: 



oc J2 l(^IV'(ro))Pe-^(WM» 



(34) 



Here, i/'(r) denotes entire microscopic field configuration 
at the imaginary-time instant r, which corresponds to a 
quantum state. The overlap |('0|'0(to))P is taken at an 
instant tq for which is specified in the Heisenberg 

picture. A fundamental property of the overlap is that it 
becomes larger when the state \ip{To)) is more "similar" 
to the state (and very sharply so if the diagonal rep- 
resentation is used to construct the path- integral) . We 
will use this property in order to simplify calculation and 
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only estimate the probability of the system being in a 
neighborhood of the state ip. Having in mind that we 
are looking for a potential long-range order, which would 
unavoidably imply existence of a static order parameter, 
we can seek probability of the system spending an ex- 
tended period of time in a neighborhood of the state ip. 
The quantity that we will consider is: 



E n 



-l(V',V'(T))e-^»'^(^)» . (35) 



The neighborhood is specified by the function 9{ip^ ip') 
that becomes smaller when the states ip and ip' become 
more different. Physically, F{ip) is the free energy asso- 
ciated with the fluctuations from the vicinity of ip. With 
a proper choice of the neighborhood function, the states 
with smaller free energy are more energetically and en- 
tropicaly favored. 

We now turn to our specific problem given by the ac- 
tion (|29|l . Recall that the saddle-point vectors ^ corre- 
spond to static spin states. Since the action is expanded 
about the saddle-point ^ in the expression (j^ . any non- 
zero values of the height fields mean moving away from 
that saddle-point. This allows us to define the neighbor- 
hood function 9 in a soft, but controlled way: 



2 T \ 

-gm X X) 



(36) 



where is a tunable parameter controlling the neigh- 
borhood size. Then, the free energy of a state is: 

e-^(«) - ^e-^'^'^^^) (37) 

X 

^ g-3[x^Cx+(x^g-l-g^x)-l-™^X^x] 
X 

Without the "mass" term rn^ the free energy would not 
depend on the saddle-points, because mere shifts of vari- 
ables in the path integral that leave the action invariant 
would switch between them. This addition to the the- 
ory does not alter its fundamental properties. The phys- 
ical bare modes are "gapped" to begin with, and the 
mass only changes the gap in the coupling matrix. On 
the other hand, the "gapless" unphysical bare modes do 
not couple to the physical degrees of freedom, and giving 
them mass is only a convenient way to handle them. Also 
note that the neighborhood function does not alter the 
nature of dispersion in the action, and thus cannot intro- 
duce unwanted effects at the neighborhood boundary. 

We can make progress in two limits: very small and 
very large coupling constant g. For g <^ I the summation 
over integer fields x in H37|) can be approximated by an 
integration: 



Vxe 



(38) 



oc exp g^^{C + m'^)-^$, . 



However, since all the saddle-point vectors are degener- 
ate eigenvectors of the coupling matrix (|33|l and have 
the same normalization (|32|) . the free energy will have 
no dependence on the saddle-points. This limit together 
with large corresponds to the calculation from the 
Reference;^ only performed to all orders of 1/m^. 

For gm^ 3> 1 it is convenient to perform the Poisson 
resummation in (|37|) . Let us introduce a vector /i with 

integer components /i^'^'' (fc — 1,2), and write: 



oc 



/ ^?Xexp[-5'(x;6-«4M^X + X^M) 

f - oo 

^exp[.g(^+— /x) (C-fm2)-i(^ + — /X 



This expression simplifies considerably due to H32|) and 
(Pt: 



-Fii) 



2^ exp 



-/x^(C + m2)-i 



(39) 



6 + 171- 



Note that we could also add an explicit "vortex core" 
term ufi^ fi to the free energy. It would soften the 
integer-valued constraints for the height fields, and yield 
a sine-Gordon theory in the large u limit. Doing this 
would be useful if we needed to discuss stability of phases, 
but for the purposes of present problem, this will prove 
to be unnecessary. The smallest eigenvalue of the matrix 
C + m? in the last equation is m? . Therefore, in the 
limit of gm^ 3> 1 (and u ^ 1), the only rapidly varying 
part of the exponent on the right-hand side is the purely 
imaginary part. We can easily understand the oscillatory 
effect that it induces as long as -C 1. First, note that 
the components of the saddle-point vectors ^ always have 
integer values: all possibilities are shown in the Table ^ 
Then, we can decompose the integer-valued components 
of the vector /x into two parts: 



(fc) 



(k) 



6fi 



(k) 



<5/.f e{0...5}. (40) 



Since the quadratic part of the exponent in (|39|l varies 
only very slowly, we can neglect fluctuations of 5fi in it. 
Similarly, since <IC 1, we can neglect fluctuations of 
M. in the oscillatory part (they approximately contribute 
a 27r X integer phase). We approximately have: 



e ^^^^ oc exp 



(Gtt) 



6 -h 



(41) 



Clearly, the oscillatory part will give rise to a destructive 
interference for every non-zero component of the saddle- 
point vector ^ (yielding factors of the order of m^, or 
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{gm?)~^ in the path-integral weight). In order to min- 
imize the free energy F{^), the saddle-point vector ^ 
should have as many zero components as possible. The 
appropriate quantity to consider is: 

np = 2nc + ria , (42) 

where n^, n;, and ric are respectively the total numbers 
of the A, B, and C type sites from the Table HTl in a 
saddle-point vector. Configurations that maximize Up 
are preferred, and entropically selected by fluctuations. 

Simple algebra can be worked out to find Up. The total 
number of Kagome lattice sites is: 

ria + niy + nc^ N , (43) 

while the total number of dimers in a least frustrated 
state is: 

1. , 2N 

-(2n, + nfc) = — , (44) 

since one dimer sits on every Kagome triangle, and the 
triangles share corners instead of bonds. Combining 
these two equations we find: 

N 2N 
"^a = y + , ^ ~3 " ■ ^ ' 

In order for all n^, rif,, and Uc to be positive and smaller 
than N, ric must be bounded between and iV/3. We 
see that Up is maximized simply when the number of C- 
sites is maximized, which means: Ua = 2iV/3, Ub = 0, 
fic — N/3. Note that the preferred configurations also 
have the maximum number of flippable spins, whose flip- 
ping costs no energy. Every A-site is a flippable spin, be- 
cause the numbers of frustrated and unfrustrated bonds 
emanating from it are equal (flipping a spin toggles bond 
energy on every emanating bond). Some untypical pre- 
ferred configurations are shown in the Fig. ITTI 

The total number of preferred configurations is macro- 
scopically large. This can be demonstrated by observ- 
ing that they map to the hard-core dimer coverings of 
the honeycomb lattice. The preferred configurations have 
only A and C type sites. Two C-sites cannot be neigh- 
bors, but their number should be maximized, so that 
every C-site can be represented by a dimer on the cor- 
responding honeycomb lattice bond, as depicted in the 
Fig. 1121 A transition graph can be found by overlap- 
ping any two honeycomb lattice dimer coverings, and it 
consists of isolated loops, the smallest having six hon- 
eycomb bonds. Therefore, the preferred configurations 
on the Kagome lattice are only locally different from one 
another, and may be transformed into one another by 
flipping six or more flippable spins (one at a time). 

E. "Disorder-by-Disorder" 

The conclusion so far is that the maximally flippable 
states may be entropicaly selected by fluctuations, to 





A.Z 


A2^ 


7 
B — 


cA 


EC.. 


±2 





±1 





E 





±2 


±1 






TABLE II: All possible local configurations of dimers (frus- 
trated bonds) at the saddle-points, and the corresponding val- 
ues of the saddle-point vector components. The site i sits at 
the center of the bowtie, and may be of type A (two vari- 
eties), B, or C. The description of the saddle-points in terms 
of dimers is given in the section 1111 CI 




(a) (b) 

FIG. 11: Two periodic (not typical) preferred dimer config- 
urations (with maximum fiippability). Every site with two 
dimers emanating from it holds a fiippable spin, and only 
every third site holds an unfiippable spin. 




FIG. 12: The maximally fiippable states map to the honey- 
comb lattice hard-core dimer coverings. The Kagome triangle 
centers form a honeycomb lattice whose bonds go through the 
Kagome sites. For every Kagome bowtie with a C element 
from the Table ^ put a dimer on the honeycomb lattice. 
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smaller or greater extent depending on the value of g. 
If we naively traced back connections between g and the 
parameters of the original spin model, we would find that, 
generally speaking, smaller g describes stronger further- 
neighbor and multiple-spin exchange processes. 

In principle, the maximally flippable states will be 
mixed together in the ground state, since it takes only lo- 
cal fluctuations to change between them. Our findings so 
far have included only the effects of "small" fluctuations 
about various saddle-points. Even though a smaller value 
of can always be chosen to expand the scope of in- 
cluded fluctuations, too small may invalidate the ap- 
proximations that made the calculations possible. Hence, 
certain "large" fluctuations are beyond reach of this for- 
malism. Exactly these fluctuations decide whether the 
maximally flippable states are evenly mixed with all other 
states or not. If they are, the ground state is chaotic and 
disordered. The only chance for a valence bond order is 
if something suppressed the fluctuations into all but the 
maximally flippable states. Then, the effective degrees of 
freedom would be only the honeycomb lattice dimer cov- 
erings, and their local dynamics would yield a plaquette 
dimer long-range order typical of dimer models on bipar- 
tite lattices. It is our goal here to determine whether 
such long-range order might be stable. 

There are several arguments that can be made in fa- 
vor of the disordered phase. The first thing to observe is 
that the minimum of the free energy is extremely widely 
and evenly distributed over a large number of disordered 
states (mappable to the honeycomb lattice dimer cov- 
erings). This means that the system does not spend 
much time fluctuating near any particular one of them, 
likely ruling out a static order parameter, and thus lattice 
symmetry breaking. Secondly, due to the degeneracy in 
the action, it is possible to make local field changes that 
correspond to flipping a single spin without paying en- 
ergy on the spatial links. Since such small fluctuations of 
the flippable spins are energetically controlled only along 
one dimension (imaginary time), there is nothing to stop 
them from proliferating. This would be true even if the 
"smooth" phase were obtained in the "height" action: 
the macroscopic degeneracy due to geometric frustration 
allows many spatially different patterns that break the 
"height" symmetry. If fact, every least frustrated state 
has a macroscopic number of flippable spins (see from 
(|45|l that Ua ^ y). Therefore, there is no mechanism 
to suppress fluctuations into any possible state. Note 
that for all but the maximally flippable states to be sup- 
pressed, the only favorable flipping processes would have 
to simultaneously involve at least six spins (dimer flip on 
a honeycomb lattice hexagon). 

Another consequence of abundant single-spin fluctua- 
tions is absence of magnetic order in the ground state. 
In our problem this also contradicts possibility of the va- 
lence bond order, since it would be accompanied by a 
net Ising moment. All maximally flippable states have 
macroscopic magnetization M = To see this, note 

in the Fig. ^2 that all flippable spins (A-type sites) must 



be aligned, since they are connected to each other either 
through one frustrated bond (dimer), or through two un- 
frustrated bonds. 

It is apparent by now that all minimally frustrated 
spin configurations are mixed into the disordered and fea- 
tureless ground state. The correlations are short-ranged 
since there is a macroscopic number of flippable spins in 
every least frustrated state, making the spins virtually 
independent. A property that distinguishes the Kagome 
from the other lattices is the formal lack of dispersion 
in the lattice field theory. We interpret this as a signal 
that the excitations are very heavy or perhaps even local- 
ized (exactly true in the small-g limit). Indeed, strictly 
short-ranged spin-spin correlations, consistent with very 
undispersive modes, have been observed in the Monte 
Carlo simulationsiS2i2S 



F. Quantum-Mechanical Interpretation 

The analysis of the lattice field theory has yielded two 
essential results that can help us sketch the ground and 
excited states of the Hamiltonian Q. They are: a) 
no symmetry is spontaneously broken, b) the excitations 
have localized character (very large effective mass). We 
use the free energy F{^) from the previous section as a 
rough indication of the probability amplitudes that differ- 
ent dimer configurations have in the ground state. The 
ground state is a smooth superposition of all possible 
configurations [ip) of frustrated bonds: 

|0)=5]a^|V). (46) 

The amplitudes of the similar states are roughly equal in 
magnitude. This is required in order for two states dif- 
ferent by a single spin flip to give a large matrix element 
{ipi\(—TS^)\'ip2) and yield a significant energy gain. How- 
ever, the amplitudes depend on flippability of the states 
IV')- The state with a larger number of flippable spins 
will have a larger probability | p . 

Due to a very localized nature of excitations, we can 
say that the physics of this model is very similar to 
the physics of completely disconnected quantum spins 
in transverse field, for which all eigenstates are known. 
Although the actual flippable spins interact, their inter- 
action seems to be largely inconsequential. This suggests 
that many good variational wavefunctions (for ground 
and excited states at F <g; J^) can be obtained by a sim- 
ple Gutzwiller's projection: take the states of the non- 
interacting Kagome spins in a transverse field and project 
them to the manifold of least frustrated states. All exci- 
tations are gapped, and the gap is ~ F. 

Finally, we recall that this disordered quantum phase 
is not topological in the original spin model. Clearly, it is 
stable against small higher order perturbations in F/J^. 
In fact, it obtains for all values of the transverse field F, 
without any intermediate phase transitionsiS2i22i2£ 
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IV. HEISENBERG MODEL WITH EASY AXIS 

In this section we analyze the XXZ model on the 
Kagome lattice, and its extensions: 



J, 



J 



SfS]) 



(47) 



with Jz,J± > 0, and J± J.^. The analysis will closely 
follow that of the transverse field Ising model in the sec- 
tion mil and rely on the notation and conventions defined 
there. Many similarities will be encountered, except that 
the formalism will be of greater complexity. One appar- 
ent difference, however, is that the total magnetization 
in z direction is a good quantum number in this model. 

As before, we begin by considering an effective theory 
that describes the physics at the energy scales well below 



= ^T.MsrSi+SrS+)Po + o{:^) . (48) 



cff 



This theory lives in the Hilbert space spanned by the least 
frustrated states of the pure Ising model, and Vq is the 
projection operator to this space. It describes dynamics 
of the flippable spin-pairs on the Kagome lattice bonds. 
In order for a pair of spins to be flippable, the spin con- 
figuration must be minimally frustrated before and after 
the pair is flipped. Note that this automatically requires 
that the two spins be antialigned fFig. ll3|l . The effective 
theory can again be expressed as a soft-core dimer model 
in the same Hilbert space as the one that described the 
TFIM model, but with more complicated dynamics. We 
will reformulate it as a U(l) gauge theory, derive a dual 
lattice field theory for it, and explore the possible phases. 
This time we will find a valence bond crystal and a spin 
liquid. 



A. U(l) gauge theory 

The only difference between the XXZ effective dimer 
model and that of the TFIM model Q is that now 
the elementary loops on which the dimers can be flipped 
enclose two dice lattice plaquettes, instead of one (since 
two spins are flipped at a time). There are four such 
processes, and they are shown in the Fig. FPil 

The U(l) gauge theory is built the same way as in the 
section IIII Al The electric field and the charge-2 bosons, 
whose fiuctuations are controlled by the Gauss' Law H12(l 



FIG. 13: Flipping a pair of aligned spins creates extra frus- 
tration on the triangle that contains them. Dimers denote the 
frustrated bonds. 



^ 

(a) 

q6 qe 

C^l ^ 

P6l'^'33^P6r P6l'^'53^P6r 

(b) 

qe qe 

1^] ^ ^ 

(c) 

qe qe 

^ 

P6 - >^P6 R P6 ly^- ->^6 R 



(d) 

FIG. 14: Four possible low energy processes that keep frus- 
tration at the minimum. Recall that the minimum of frustra- 
tion is achieved if there is one dimer emanating from every 
3-coordinated dice site, and an even number of dimers em- 
anating from every 6-coordinated site. The dashed arrows 
show between which two 6-coordinated sites a pair of dimers 
is exchanged (charge-2 boson hopping in the U(l) gauge the- 
ory). These processes preserve the global Ising magnetization: 
the bond between the two plaquettes is always unfrustrated. 



and the potential energy (jSJ, represent the low energy 
degrees of freedom. The new form of the kinetic energy 
can be easily obtained by comparing the two-plaquette 
processes in the Fig. with the single-plaquette pro- 
cesses in the Fig. El The two single-plaquette processes 
consistent with the low energy physics can be combined 
in four different ways to give the allowed two-plaquette 
processes. In combining them, the middle bond {qeqs) 
is fiipped twice, so that there is no net change on it. 
This gives us the operators in the Table IIIII The ar- 
gument of each cosine is the sum of two corresponding 
single-plaquette circulations and boson hopping(s) from 
the expressions © and (|10|l . but multiplied by the fac- 
tors of r](^pq^ and e(pg) in such a way that the contribution 
of the central bond (qeQa) is properly canceled out. Then, 
if we label the four processes by a = a,b,c,d, the U(l) 
effective Hamiltonian is; 



(pq) ' 



(a) 
kin 



(49) 
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J± Yl cos(/i 



(<36 93)' 



J± J2 cos (ft'"' 



-E- 



(gels)' 



(legs)' 



^1" . 

(96 93) 



(96 93) 



(96 93) 



■ ^(9693) E '-^P9 + ^(9693) E ■^P9 

P9 P9 

03 / c% 

■ ^(9693) E -^P9 + I (¥^96 ~ VP6H ) '?P6fi 96 + Yl ^ (pq) -^P 

pq \ pq 

■ ( iVpeL ~ fqs}''lqiiP6L + E £(P9)'^P9 ) + ^(9693) E ■^P9 
V P9 / P9 

0^ 



I (VP6L ~ fq&^'HqaPSL + E £(p9)'^P9 ] + I ('y^96 ~ '/'p6h)^P6H96 + E ^(p9)-^P9 1 

\ P9 / V pq J 



TABLE III: Kinetic energy operators corresponding to the processes in the Fig. 1141 Notation for the sites is defined in the 
figure: (qeqs) is the bond shared between the two dice plaquettes, while Pql and p^R are the bottom 6-coordinated sites on the 
left and right plaquette respectively. Only the bonds with arrows are included in the sums. 



The U ^ oo limit is an exact rewriting of the effective 
dimcr model, while a finite U introduces various new dy- 
namical processes, defined on larger dimer loops (spin 
clusters), but consistent with the global spin- flip symme- 
try of the XXZ model. 



In terms of the fields -^[pg^: the magnetic field and the 
particle current (labeled by the Kagome sites) are: 



pq 



}_^Jpqj^(b) 
2 (P9> 



l±Jpq_j^{c) 

9 (pg) 



B. Lattice Field Theory 

We proceed by writing the path-integral for the U(l) 
Hamiltonian, respecting the constraints given by the 
Gauss' Law (|12f) . The action will contain the Berry's 
phase (|13|l and the potential energy term p4|) as be- 
fore. However, this time the Villain's approximation 
will give rise to four integer-valued massive fields ^[p^^ 
{a = a, b, c, d) that live on the dice lattice bonds, and 



couple to the arguments of cosines /i^^^^ 

nn 



from the Table 



'kin 
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(pq) V (pq) 



An additional Berry's phase ~ ittK''^^^ appears because 
the coupling J_l is positive. It is again possible to define 
the magnetic field Bi and the particle current jpgqg- The 
magnetic field couples to the plain plaquette curl of the 
vector potential, while the current couples to the boson 
hopping: 
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We sketch the derivation of these equations in the Fig. IT^ 
and its caption. The magnetic field and the particle cur- 
rent are integer- valued fields, and the action is seemingly 
reduced to the form that pertains to the TFIM model. 
When the vector potential Apq and the particle phase Lp^^ 
are integrated-out, the same current conservation 
and Maxwell's (|18() equations are recovered. 

We can again introduce the height fields Xi Xi as 
in 1(201). For this purpose, it would be convenient to ex- 
press the fields as time derivatives. Since the path- 
integral has a closed boundary condition in imaginary 
time, and ^/pi) t ^^ed not be zero, we can generally 



write: 
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where the time dependencies have been explicitly shown 
(/3 ^ oo is the "last" moment of time). Ignoring for a 
moment the magnetic field and current terms from H5U|I . 
the kinetic part of the action becomes: 



'S'kin 



E 

{pq),a 



■/3-1 



E»(^-'«)' 



(53) 



.T = 



(pq)fi 



''{pq),l3 



K 



(a) 

{pq) 



ittK 



(a) 

{pq) 



E 

{pq),a 



.r=0 



(a.) _ (a) ^ 
«'^"^{pg),0 '^{pq),P. 



15 




(a) 



(b) 



(c) 



(d) 



FIG. 15: Explanation for the derivation of equations (I5H : 
The goal is to associate Bi to the curl of Apq, and ji to the 
curl of e^pq^Apq. For each of the four plaquette pairs in the 
figures (a) through (d), a counter-clockwise circulation (of 
Apq or £{pq)Apq) IS takeu on the central plaquette, labeled 
by its dual Kagome site i. The fields K^^^^^ on the empha- 
sized bonds will be coupled to this circulation, and we simply 
collect those coupled to the curl of Apq into Bi, and those 
coupled to the curl of e^pq)Apq into ji. However, there is a 
small complication. Let us call a plaquette "left" or "right" 
according to its position when the plaquette pair is rotated 
to point "upward" like in the Fig. 1141 Then, for the pairs 
(a) and (c), the plaquette i is "left", while for the pairs (b) 
and (d), i is the "right" plaquette. The fields K^^'^^^ and ^^^p^^ 
couple to circulations of the different quantities on the "left" 
and "right" plaquettes (see Table IllH . Therefore, the "left" 
and "right" must be distinguished. Notice that the oriented 
emphasized bond in the figures (a) and (c) points from the 
3-coordinated site toward the 6-coordinated site {ripq = — 1), 
and the opposite in the figures (b) and (d) {rjpq = +1). This 
can be used to determine when the circulation is made on the 
"left" or on the "right" plaquette, and this is the origin of the 
r)pq terms in USUI . 



would be suppressed. The macroscopic physics of such 
triplet bonds must be very similar to that of the singlet 
bonds, because in both cases every spin can be paired 
with only one of its neighbors. Even the same higher 
order processes, such as the valence-bond movements on 
the closed loops, would be preferred by either sign of J±. 

Passing completely to the Kagome lattice notation, 
and neglecting the Berry's phase, the action of the final 
field theory becomes: 



(54) 



where the fluctuations of the height fields are constrained 
by the equations (|51|l in the dual form: 
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(a) 1 - VtJ Jb) , 1 + Vtj Jc) 
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(55) 



C. Important Properties 



The field K^^^^ was integrated out in the last line, and an 
emerging additive constant was discarded. The bound- 
ary conditions in imaginary time now appear open, and 
the extra Berry's phase, due to the positive J±, appears 
only at the boundary. In fact, the Berry's phase is sensi- 
tive only to the parity of the integer-valued fields at the 
boundary. Whenever the fluctuations render this par- 
ity short-range correlated along time, one may expect 
that the Berry's phase will not affect the macroscopic 
properties of the theory. This will certainly happen in 
any disordered phase. However, it can also happen in a 
"smooth" phase that describes a plaquette valence-bond 
order: the macroscopic degeneracy created by geometric 
frustration allows many locally different "smooth" states, 
and the small fluctuations between them are extremely 
abundant, especially on the corner-sharing lattices such 
as the Kagome. 

Therefore, we will assume in the following that the 
Berry's phase shapes only certain local properties of fluc- 
tuations, and neglect it for the purposes of discussing the 
possible phases of the theory. The affected local proper- 
ties can be revealed from a microscopic point of view. 
The positive value of the coupling J± in the XXZ model 
(|47|l prefers the spin singlet formation on the Kagome 
bonds: \ H) — \ |t)- If J± were negative, the symmetric 
triplets would be favored instead: | Ti) + I it)- Note that 
at least due to the strong Ising antiferromagnetic interac- 
tion, the "true ferromagnetic" nature of the negative J± 



The action H54() resembles very much that of the TFIM 
model (|27|l . However, it could give rise to a very differ- 
ent physics. Certain kind of fluctuations are forbidden by 
the action H54II , and the remaining ones might be able to 
entropicaly select some ordered state. The forbidden fluc- 
tuations are those that change the total magnetization of 
the Kagome Ising antiferromagnet. The mechanism for 
this is provided by the constraints on the spatial con- 
figurations of Xi and Xi, which emerge from (|55ll . The 
natural degrees of freedom in this field theory are j 
and they live on the Kagome bonds, reflecting the nature 
of the XXZ perturbation. 

Therefore, we will adapt the analysis of the TFIM 
model to these new degrees of freedom, and write the 
action (|54|l in the matrix form. First, we note that 
the action is minimized by the same height field config- 
urations Xi ^nd Xi as before. By shifting variables and 
expanding about a particular saddle-point, the potential 
part of the action becomes: 

Spot = 5 E E ~ + ^iv) -^3)+ ^^j) ' (56) 

where has been defined in the section IIII CI Then, 
we apply the resummation formula: 

E E ^(y) = E(^(*J')"* + (57) 

i j£i {i]) 
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(Cpot/t) '"I = £{ij) Uxi - E (Xfc + e<ife) Afc)) + £<ij> Uxi - E (Xfc + eofc>Afe)) |S"?> = £(,j) E Cifc + e(ij> E ^jk 
^ ^ kei ' ^ kej ' kei kei 

(CpotK)*^]^ = [(^4Ai - Y.{^k+ e(^k)Xk)) +e<ij> (4xi - E (Xfc + eofc) ^fc))] , if^j) = [E e<ifc>?»fc + e<ij> E Ofe] 
(CpotK)^^' = Li^) ('4xi - E (Xfe + + (4Aj - E (-^fc + eofc>Xfe)) ^S-]) = \ E Cife + E eOfc>?ifc 

{Cpotit)[f = ('4A, - E (Afe + e(»fe>Xfe)) + (4Aj - E (^fe + eo-fe>Xfc)) = E e(!fe>^»fe + E £{jk)^jk 

TABLE IV: Action of the spatial (potential) part of the coupling matrix on the field vectors (left), and the components of the 
saddle-point vectors (right). For the components of type (b) and (c), the sites i and j must be chosen in such a way that the 
bond orientation is from i toward j. 



to the expressions (|55|l . and substitute the result in (|30|l . 
This gives us the matrix form of the action H54|) : 



(58) 



Components of the vector k are the k\"^^ fields. The 
structure of the coupling matrix in terms of the natural 
degrees of freedom, and the saddle-point vectors are given 
in the Table |IV| 

Now, we repeat the analysis from the sections IIII CI 
and IIII Dl in order to find the effect of fluctuations. The 
crucial pieces of information are how the saddle-point 
vectors are normalized, and how the coupling matrix acts 
on them: 

• all saddle-point vectors ^ have the same norm: 

= const. ; (59) 

• all saddle-point vectors ^ are degenerate eigenvec- 
tors of the coupling matrix C: 



(60) 



In full analogy to the TFIM case, the coupling matrix C 
is completely dispersionless. There are 24 bare modes per 
unit-cell of the Kagome lattice (unit-cell has six bonds), 
and all of them are localized. Only four of them have a 
non-zero eigenvalue at zero frequency (equal to 36), while 
the other 20 are "gapless" and unphysical fluctuations 
(due to the redundancy of representation). 



D. Effect of fluctuations 

In the quest for potential order-by-disorder, we proceed 
in exactly the same fashion as before. The XXZ model 
only brings a new complication: the total Ising magneti- 
zation is conserved. Every value of total magnetization 
defines a separate sector of states, and the fluctuations in 
the lattice field theory (|54() cannot mix the states from 
different sectors. In principle, the entropical effects of 
fluctuations should be investigated for every sector sep- 
arately. However, we only need to focus to the sector of 



zero Ising magnetization, since the XXZ coupling clearly 
favors it. 

We introduce the free energy F{^) of fluctuations 
about a saddle-point ^, and search for the saddle-points 
that minimize it. For small g the same situation occurs 
as in the TFIM case (see equation pSf) ): there is no 
entropical selection of the saddle-points. Thus, a disor- 
dered ground state is obtained, which in the XXZ case in 
fact has non-trivial topology, as will be argued at the end 
of the section IIV b'\ New interesting things happen for 
large g. The free energy F{^) in the large g limit, after 
the Poisson resummation, is given by the approximate 
expression analogous to (|41() : 



^exp 



(36^)^ 



36 -|- vn? 



where the integer Poisson fields P'^!^^^ {a — a . . .d), form- 
ing the vector /i,, have been decomposed as: 



(61) 



(q) 



.35}. (62) 



The free energy is minimized when the saddle-point vec- 
tor ^, given in the Table |lVl has the maximum number 
of zero components. All other possible values of the com- 
ponents ^ are integers and factors of 36, so that only the 
zero components avoid destructive interference in H61|) . 

In order to discover which saddle-point dimer cover- 
ings are preferred and minimize the free energy, we need 
to characterize them in terms of the local dimer configu- 
rations at the Kagome sites, bonds and triangles. The 
Kagome sites have already been characterized by the 
number of dimers emanating from them, in the Table 
im All non-equivalent dimer arrangements in the neigh- 
borhood of a Kagome bond are systematically shown in 
the Table together with the corresponding numbers 
of the zero components of the saddle-point vector. Fi- 
nally, the triangles can be characterized by the types of 
sites at their corners, and all possibilities are given in 
the Table IVII For every allowed type of triangle one 
can find three situations in the Table corresponding 
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B, C 




aYb 




i^B 


4 





2 





2 





2 



TABLE V: All possible non-equivalent configurations of frus- 
trated bonds in the neighborhood of a Kagome bond. The 
types of sites are labeled according to the scheme from the 
Table [H] (for the bottom site on the central triangle there 
are always two options). This table shows the number of 
zero-valued components , a = a ... d of the saddle-point 
vector, for the horizontal bond {ij} on the central triangle. 



to the three bonds on the triangle (one of which is frus- 
trated), and collect the total number of the saddle-point 
zero components that such a triangle would contribute. 
Adding contributions of all triangles, that is all bonds, 
we obtain the following "scoring" number that should be 
maximized: 

Tip = Subbb + '2nabb + "^riaab + ^riaac + Z^abc + ^nbbc ■ (63) 

The quantities Ubbb ■ ■ ■ n-bbc denote the total numbers of 
various kinds of triangles in a given saddle-point dimer 
configuration. At this stage, we have to investigate possi- 
ble relationships between these numbers. The first thing 
to note is that the total number of Kagome triangles is: 



nbbb + nabb + naab + naac + nabc + nbbc 



2N 
~3~ 



(64) 



Then, using the Table IVII we can count the total num- 
bers of A, B, and C sites: 

"a = ^C^riaab + '^Uaac + "abb + nabc) , (65) 
m = -[Subbb + '2nabb + '2nbbc + naab + nabc) , 
nc = -^{naac + nabc + nbbc) ■ 

By combining these equations with l|45(l . one finds that 
Ua, Ub, nc, Uabb, and riabc can be expressed in terms of 
the independent variables Ubbb, naac, nbbc, and Uaab- The 
"scoring" number n'^ can now be simplified using the 
identity (|64(l . The quantity that has to be maximized 
is: 



Up = iUbbb + '2na 



nbbc 



(66) 



and the variables appearing in it are independent, al- 
though subject to inequalities ^ ria ^ N , . . . , ^ 
nbbc s; 2N/3. 

Finding the absolute maximum of rip is a well posed 
problem of linear programming. It can be shown that it is 
obtained when the number of C-type sites is maximized, 
which gives rise to the very same states preferred by the 
fluctuations of the TFIM model (see Fig. Ill|l . However, 
these are not the saddle-points that we are looking for: 




Aa a JX 3ZX 



BBB 



ABB AAB AAC ABC 



BBC 



TABLE VI: Characterization of all possible kinds of Kagome 
triangles. The site types A,B, and C are defined in the Table 
im Note that two C-sites can never be neighbors, and that 
three A-sites cannot sit on the same triangle. The dimers 
emphasize which bond must be frustrated in a given situation 
(there are multiple choices only for the case BBB). 



they have a macroscopic magnetization. Besides, instead 
of being maximally flippable, as in the case of the TFIM 
model, they are now minimally flippable; in fact, they 
have no flippable spin-pairs at all. If we want to find 
the configurations with zero magnetization that maxi- 
mize Hp, we must explore a path different from having a 
large number of C sites. This is, in principle, a difficult 
problem, and an exact analytical solution is not avail- 
able at this time. Instead, we guess that Ubbb should be 
made as large as possible. This yields the configurations 
without any C sites, and with a large number of B sites 
(see (gSl)- The best choice is: Ua = N/3, rib = 2N/3, 
nbbb — naab = N/3. It is possible, though relatively 
complicated, to demonstrate that the configurations with 
these parameters break the lattice symmetries in a unique 
stripe-like fashion shown in the Fig. ^| It turns out 
that such states are not magnetized at all, and that they 
have a large number of fiippable spin-pairs. Their scor- 
ing number (rip = N) is significantly larger than that of a 
typical unmagnetized state (a fraction of N) . Therefore, 
they are excellent candidates for the preferred configura- 
tions. No better configurations were found when every 
least frustrated state with zero magnetization was ex- 
plicitly examined using the computer (the sample had 24 
sites and closed boundary conditions). 

In conclusion, for large g the preferred configurations 
of frustrated bonds that minimize the free energy break 
the lattice symmetries in the stripe-like way, as shown in 
the Fig. El This indicates that a valence-bond ordered 
phase could be realized in the XXZ model when dynam- 
ics is dominated by the short-ranged spin-pair flips (the 
larger g, the weaker further-neighbor and multiple-spin 
exchange) . It is now necessary to verify stability of such 
an ordered phase. 



E. Stability of the Valence-Bond Order 

A usual way to determine whether fluctuations ulti- 
mately destroy the long-range order involves rcnormal- 
ization group (RG). The lattice field theory H54|l of the 
XXZ model resembles an integer-valued height model, 
and one might naively hope that the RG arguments could 
be applicable to it. In a standard and simple integer- 
valued height model (on the square lattice, for example) 
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FIG. 16: The preferred saddle-point configuration witfi zero 
magnetization, selected by the XXZ fluctuations. The lattice 
symmetries are "broken" in the stripe-like fashion. There are 
only type A and B sites in this state, and the number of 
BBB triangles is large (they sit between the straight chains of 
dimers). The total magnetization is zero: every two spins con- 
nected by a vacant bond are antialigned, so that the straight 
chains alternate in magnetization, as well as the dimers in the 
middle along the chains. 

one first softens the integer constraints for the height 
fields by writing a sine-Gordon theory. Then, one checks 
how the sine-Gordon coupling flows under RG, starting 
from various parameter values in the theory. If it flows 
toward zero, then the integer constraints are irrelevant 
at the macroscopic scales, and the height model may be 
found in the "rough" disordered phase. Alternatively, the 
flow can be toward infinity, in which case the "smooth" 
long-range ordered phase is realized. In the context of 
frustrated magnetism, the appropriate height model typ- 
ically comes with a background in our case), so that 
the "smooth" phase also breaks the lattice symmetries»2i 
Therefore, let us write a sine-Gordon theory for the 
XXZ model, based on the action if^ : 

The K/".^, fields are now real-valued, and their deviation 
from integers is penalized by the sine-Gordon term, es- 
pecially in the large 7 limit. It is convenient to shift 
the variables k by C^$, = ^/36, and remove the linear 
terms: 

T,{ij),a \ / 

Unfortunately, it is not possible to directly apply the 
RG treatment to this theory. The bare modes (modes of 
the coupling matrix C) are not only dispersionless, but 
some of them appear gapless as well. The "gapless" bare 
modes are redundancy of representation, but they still 
pose a technical difficulty. It is through the sine-Gordon 
coupling that they at least acquire dispersion. The sine- 
Gordon term mixes the bare modes when they describe 



non-integer fluctuations of the height fields. Let us rela- 
bel the fields as where r is a vector specifying 
a Kagome lattice unit-cell, and n G {1,2... 24} is an in- 
dex specifying the bond {ij) (one of six) and the flavor a 
(one of four) within the unit-cell. Then, we can express 
the fields as the linear combinations of the 24 local bare 
modes ^n,r- 

24 

Kn,r = Wnm,Ar^rn,r+Ar ■ (67) 

m— 1 At* 

One can formally integrate out the four physical (mas- 
sive) modes: $2i,r, *i>22,r, *i>23,r, <I*24,r, and obtain the ef- 
fective theory as a perturbative expansion in 7. The ef- 
fective theory is a complicated expression involving co- 
sine terms whose arguments are linear combinations of 
the remaining twenty modes: 

20 

^off = sEE^^-"^"^-)" (68) 

n— 1 r 

(20 
m— 1 Ar 
-2^^Cn,r^ + 0(7') • 

The redundancy of representation survives in the effec- 
tive theory through periodicity of the cosines. However, 
the redundancy is easily removed by treating the mode 
amplitudes in the effective theory as angles: the physi- 
cal degrees of freedom (which have been integrated-out) 
enter the effective theory precisely trough the residual 
[0, 27r) amplitudes of the remaining modes. 

Let us for a moment expand the cosine terms to the 
quadratic order, and obtain a Gaussian theory for the 
twenty modes. In absence of the saddle-point background 
^ = 0, such a theory would be trivial: 

This means that the effective theory (|68ll retains the 
nature of the original sine-Gordon model: small fluctu- 
ations of the bare modes are gapped. However, we are 
concerned only with the vicinity of the candidate valence- 
bond ordered state (Fig. [TB|) . For the saddle-point vector 
^ that describes this stripe pattern the effective band 
structure in the Gaussian approximation is shown in the 
Fig. 1171 Apart from having frequency dependence, the 
lowest lying mode is dispersive in the direction along the 
stripes, but not in the perpendicular direction. Never- 
theless, full spatial dispersion can be expected if one goes 
beyond the Gaussian approximation, because the lowest 
lying mode is actually coupled to some higher modes that 
are dispersive in the perpendicular direction. The results 
of this approximation are only good for arguing that dis- 
persion emerges. There is no simple way of telling what 
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FIG. 17: Effective band structure of the sine-Gordon tfieory 
1681 1 for the candidate valence-bond order in Fig. 1161 (zero 
frequency, quadratic approximation). The wavevector q is 
taken (a) along the stripes, (b) perpendicular to the stripes. 
The vertical scale is given in arbitrary units, but proportional 
to 7. The Kagome lattice unit-cell had to be doubled, so 
that 40 modes are shown; only the dispersionless branches 
are degenerate. Note that the lowest lying mode is dispersive 
in only one spatial direction, but also that some higher modes 
have dispersion in the other direction. 



changes beyond the quadratic approximation. Therefore, 
the correct way in which the lattice symmetries may be 
eventually broken may not be possible to guess from this 
information. 

We can now write the effective theory in a form 
that manifestly separates the dispersive and sine-Gordon 
parts. If we partially expand the cosines from l|68|l in 
the following way (supressing the n and r indices of S^n,r 

and Xn,r = X/ Wnm,l^r^m,r+Ar)'- 



^C„cos(^27rx- ^) = ^C„|27rxsin(^^) 4- 

acos(27ra;)+ [cos(H) - a] (l - I?^) } + , 
then for a proper choice of the constant a the effective 



theory becomes: 
5eff « 5*^Ccfi(0* - lih^ii)'^ + ^^Hi)) (69) 
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Cn COS I 27r ^ ^ WnmAr^ra,r+Ar 



m— 1 Ar 



where the coupling matrix Ccff(^) collects all space-time 
dispersion (brought up by the non-trivial ^ ^ 0), and has 
no gap at q = 0. The remaining cosines open up a gap 
for small fluctuations, thereby justifying the quadratic 
expansion that took place. This is an effective sine- 
Gordon theory. As a matter of principle, the RG treat- 
ment is now applicable. Even if the dispersion were ulti- 
mately created only in one spatial dimension (along the 
stripes), combined with the dispersion in time it would 
give a "smooth" phase for sufficiently large g and 7. If 
the full (2 -I- 1)D dispersion were obtained, then only the 
"smooth" phase would exist, since the sine-Gordon cou- 
pling would always flow toward infinity. In any event, 
existence of the "smooth" phase means that the valence- 
bond long-range order for large g is stable. For small 
5, however, the "smooth" phase is disordered since the 
fluctuations cannot select an ordered state from the de- 
generate manifold. 

This concludes our discussion of the lattice field theory. 
In the following, we take a completely different point of 
view, and provide a more physical picture of the discov- 
ered XXZ phases. 



F. Nature of the XXZ Phases 

In this section we use some simple physical arguments 
and show that the short-range valence-bond picture ap- 
plies extremely naturally to the XXZ models on the 
corner-sharing lattices. This will allow us to identify a 
physical "order parameter" for the valence-bond ordered 
phase of the Kagome XXZ model, and construct qual- 
itatively good variational wavefunctions. Also, we will 
argue that the disordered phase of the Kagome (and any 
other) XXZ model has non-trivial topology. 

Let us seek variational ground states of the Hamilto- 
nian (|2Jl that are described in terms of the singlet bonds. 
The energy minimum requirements shaped by the Jz and 
terms can be met by following these criteria: 

• number of frustrated bonds is minimized, 

• total Ising magnetization is zero, 

• number of flippable spin-pairs is maximized. 

First, we explore the circumstances in which a pair of 
spins on a Kagome bond is flippable. The XXZ pertur- 
bation ~ in ^ can flip a pair of antialigned spins, 
but one must make sure that both the initial and final 
states are minimally frustrated. The Fig. ^| shows a 
flippable pair of spins. Regardless of whether the two 
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FIG. 18: A flippable pair of spins sits on the horizontal bond 
of the central triangle. The pairs of anticorrelated spins are 
emphasized. 



antialigned spins on the central horizontal bond are in 
one or the other state, every triangle has exactly one 
frustrated bond, which is a condition for minimum frus- 
tration. This requires that the two opposite bonds on 
the neighboring triangles hold a pair of antialigned spins 
each. Clearly, energy will be gained by allowing the flip- 
pable pair of spins to resonate between the two possible 
states and form a singlet bond. If the other two pairs of 
antialigned spins were also flippable, more energy could 
be gained by turning them into the singlet bonds too. 
The attempt to create as many flippable pairs as possi- 
ble leads to the hard-core dimer coverings of the Kagome 
lattice, where every dimer represents a singlet bond (in 
contrast to the earlier representation, when the dimer 
was a frustrated bond). The corner-sharing structure of 
the Kagome lattice makes this picture extremely natu- 
ral by allowing the singlet pairs to be close-packed. The 
easy-axis anisotropy clearly shapes them as good degrees 
of freedom. 

However, the hard-core dimer coverings of the Kagome 
lattice are not quite acceptable. It is known that they un- 
avoidably have a fixed number of so called defect triangles 
that hold no dimers on their bondsiSS*^ This means that 
there would be macroscopically many triangles with all 
three bonds occasionally frustrated, and the first crite- 
rion would be violated. An example of a defect triangle 
is shown shaded in the Fig. I19r a). In order to make 
sure that two of its bonds are always unfrustrated, two 
of its spins must be always anticorrelated. One varia- 
tional mechanism that achieves this is the following. Let 
us fix the relative antialigned states of two spins on the 
defect triangle, and denote those two spins by a dimer, 
like in the Fig. Iiyt b). This allows at least the dimer on 
the top triangle to resonate as a singlet bond, but there 
are now four spins connected by the dimers, thus anti- 
ferromagnetically correlated. Certainly, those four spins 
could also resonate together while keeping their correla- 
tions, but the energy gain would be much smaller than 
that brought by a singlet bond (this is a higher order 
process). Generally, two singlet bonds are lost for every 
defect triangle. The dimers no longer represent only the 
singlet bonds, but any pair of anticorrelated spins. 

This situation can be improved. It is possible to ar- 
range the defect triangles close to each other in such a 
way that they share the singlet bonds that are going 
to be lost. Consider a so called perfect hexagon in the 
Fig. l20r a). It holds three dimers on its bonds, and there- 
fore has three defect triangles around it in a hard-core 
dimer covering. By putting three extra dimers on the 




(a) (b) 

FIG. 19: (a) A defect triangle (shaded) holds no dimers on its 
bonds. As the neighboring singlet bonds fluctuate indepen- 
dently, all three spins on it are occasionally aligned, making 
all three bonds frustrated, (b) A dimer is placed on the defect 
triangle, relaxing its frustration at all times, but simultane- 
ously correlating antiferromagnetically a group of four spins. 




(a) (b) 

FIG. 20: (a) A perfect hexagon holds three dimers, and has 
three defect triangles (shaded) around it in a hard-core dimer 
covering, (b) Covering all bonds of the perfect hexagon by 
dimers removes the defect triangles, and correlates antifer- 
romagnetically all spins on the hexagon. Only three singlet 
bonds are lost per three defect triangles. 



hexagon, the defect triangles are removed, and the six 
spins on the hexagon are forced to be antiferromagnet- 
ically correlated. Only three singlet bonds are lost per 
three defect triangles. This is clearly energetically favor- 
able, and consequently the singlet bonds will arrange in 
a way that maximizes the number of perfect hexagons. 
Every perfect hexagon can then gain additional energy 
by correlated fluctuations of its six spins. Furthermore, 
the groups of singlet bonds may be able to collectively 
resonate on the closed resonant loops. There are many 
variational states that maximize the number of perfect 
hexagons, and we show some characteristic examples in 
the Fig. El 

This variational picture seems to apply very well when- 
ever the physics of the Kagome spin models is describable 
by the short-range valence bonds. We now apply it to the 
ordered phase of the XXZ model. The lattice field theory 
in the preceding sections was able to establish stability 
of a valence-bond order without full precision in deter- 
mining how the lattice symmetries should be ultimately 
broken. It only produced information on which partic- 
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(a) (b) 

FIG. 21: Characteristic variational ground-states: (a) stripe 
pattern, (b) honeycomb pattern. The number of emphasized 
perfect hexagons is maximized (1/6 of all hexagons). In the 
honeycomb pattern there are star-shaped resonant loops of 
singlet bonds, one sitting inside every honeycomb super-cell. 
There are two possible singlet bond arrangements on every 
loop, and more energy can be gained by resonant fluctuations 
between them. 



ular microstate is most frequently visited by the system 
(pattern of frustrated bonds in Fig. ll(j|l : the other nearby 
microstates are not energetically suppressed, so that they 
are visited extremely often as well. We can now combine 
this information with the variational states. The assump- 
tion is that the correct ordering pattern must be "syn- 
chronized" with the state most frequently visited by the 
system. The Fig. compares the entropically preferred 
stripe-like configuration of frustrated bonds with the only 
two compatible variational states. The scenarios (a) and 

(b) involve the stripe pattern of perfect hexagons. Intu- 
itively, the case (a) does not seem likely, because the two 
stripe orientations of the compared states are different. 
In the case (b), the stripe orientations of the two com- 
pared states are the same, but the overlap period in the 
direction perpendicular to the stripes is relatively large. 
Perhaps the best match is accomplished in the scenario 

(c) . The honeycomb pattern of perfect hexagons is also 
the most symmetric among all variational states, and the 
set of its symmetries is the closest to that of the Kagome 
lattice. Unfortunately, at this stage there is no unam- 
biguous way of telling which match is the best, but we can 
also note that exactly this honeycomb pattern emerged 
from the analysis of the isotropic Heisenberg mode^ as 
a likely ground-state. Based on all these arguments, we 
can propose the honeycomb variational state as the most 
likely qualitative description of the ordered ground state 
of the XXZ model. 

Let us now turn to the disordered phase. Our goal 
is to show that conservation of the total Ising spin has 




(c) 

FIG. 22: Overlap between the preferred configuration of frus- 
trated bonds and the variational states. Dimers represent the 
frustrated bonds, and the emphasized hexagons are perfect, 
(a) and (b) demonstrate two possible ways of overlapping the 
stripe variational states, and (c) demonstrates the overlap 
with the honeycomb variational state. The repeating unit- 
cell of the overlap has 18 sites in (a), and 36 sites in (b) and 
(c). 



profound consequences for the topological properties of 
disordered phases. Consider an arbitrary Hamiltonian 
that is invariant under the global spin-flip. It can be 
always expressed as a sum of local Hermitian operators 
that flip an equal number of "up" and "down" pointing 
spins. The eigenstates of all such operators either have 
any particular spin on the lattice fixed, or involved in a 
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group of an even number of coherently fluctuating spins 
(for example, a group of two neighboring spins is a sin- 
glet or a symmetric triplet valence bond, appropriate for 
the pure XXZ dynamics). These eigenstates are some- 
how eventually superposed to give the ground state of 
the Hamiltonian. If all the spins fluctuate in the ground 
state so that there is no average magnetization on any 
site, then the ground state is a superposition of only the 
"valence-group" states in which every spin belongs to 
a finite even-sized group of coherently fluctuating spins. 
This is a generalization of the valence-bond states (whose 
superpositions would yield the singlet ground states). 

It is now possible to define topological sectors of these 
"valence-group" states. Choose an arbitrary pairing of 
spins within every group: every spin must be paired 
with one other spin (need not be a neighbor). Visual- 
ize the pairings by strings on the lattice that connect 
the paired spins (overlaps and shapes of the strings do 
not matter). A transition graph between any two string 
configurations can be constructed by overlapping them, 
in analogy to the hard-core dimer coverings. Then, any 
two "valence group" states from the superposition that 
forms the ground state will have the transition graph 
composed of finite closed loops (as long as the Hamilto- 
nian has only local dynamics) . If now the lattice is placed 
on a torus, there will be two topologically non-equivalent 
closed paths that go around the torus and intersect the 
bonds of the lattice. The topological sector of a string 
configuration is determined by the parities of the number 
of strings that each of these paths intersect. Two string 
configurations will belong to different topological sectors 
only if any of the paths intersects their transition graph 
an odd number times. Clearly, this can never happen if 
the transition graphs always consist of the finite closed 
loops: the ground state has a definite topology. 

Therefore, any disordered state of the XXZ model is 
automatically a spin liquid, with four degenerate ground 
states on a tours. A characteristic feature of the Kagome 
lattice (and other corner-sharing lattices) is a manifestly 
weak dispersion in the far limit for which the spin liq- 
uid obtains. This indicates that the correlations in the 
spin liquid (away from the critical point) must be strictly 
short-ranged, virtually vanishing beyond a few lattice 
constants. 



V. DISCUSSION 

We have explored two kinds of the Kagome lattice 
quantum Ising antiferromagnets. The first kind was en- 
dowed by spin dynamics that did not conserve the total 
Ising spin, and was represented by the transverse field 
Ising model (TFIM). The second kind conserved the total 
spin, and its simplest form was given by the XXZ model. 
Both TFIM and XXZ models contain only the shortest- 
range dynamical processes consistent with the required 
symmetries, acting as small perturbations to the pure 
Ising model. The considered extensions include further- 



neighbor and multiple-spin exchange dynamics, and thus 
they may reflect physics of the TFIM and XXZ mod- 
els with stronger dynamical energy scale in comparison 
to the Ising interaction. The quantum phases found in 
these models are summarized in the Table IVIII 

The disordered phase of the TFIM and related models 
was found to have trivial topology. Consistently, the ta- 
ble indirectly suggests that the same phase should be re- 
alized for all values of the transverse field. Our approach 
allowed us to gain some information about the proba- 
bility amplitudes that various spin configurations have 
in the ground state superposition. Together with this 
information, the finding that excitations appear heavy 
or localized even for weak transverse fields suggested 
the following variational wavefunctions for many states: 
the eigenstates of the decoupled spins in transverse field 
should be projected to the manifold of minimized Ising 
frustration. It is evident that the corner-sharing struc- 
ture of the Kagome lattice is responsible for the very 
weak dispersion or perhaps even localized nature of fluc- 
tuations. 

Much richer physics is found when the total Ising spin 
is conserved. The XXZ and related models give rise to at 
least two non-trivial phases. The calculations indicated 
that the valence-bond order was most likely to be found 
for short-ranged and small dynamical perturbations, such 
as the one in the Heisenberg model with strong easy-axis 
anisotropy (simple XXZ model). Furthermore, a combi- 
nation of arguments led to essentially two most probable 
ordered states, the stripe and honeycomb shaped pat- 
terns (Fig. I2J • While no good arguments to rule one of 
them were provided, we suspect that the more symmetric 
honeycomb pattern is realized in most typical situations 
(no specially favored dynamical processes). This result is 
of potentially great importance, because the same type 
of lattice symmetry breaking has been proposed to occur 
in the ideal isotropic Heisenberg model, and accounted 
for the seemingly gapless band of singlet excitations ob- 
served in numerics. The physics of the ideal Heisenberg 
model is still largely mysterious: its ground state could 
be a spin liquid. Indeed, our calculations indicate that as 
the complexity of dynamical processes increases, which is 
similar to what happens when the amount of easy-axis 
anisotropy is reduced in the XXZ model, a phase transi- 
tion into the spin liquid must occur. If our calculations 
could indeed qualitatively describe the anisotropy reduc- 
tion, a question would arise whether the phase transition 
happens before or after the full isotropy is reached. In 



dominant dynamical 
processes: 


simple 
short-ranged 


multiple-spin and 
ring exchange. . . 


do not conserve ^ SI 

i 


disordered 


disordered 


conserve ^ 5f 

i 


valence-bond 
crystal 


spin liquid 



TABLE VII: Quantum phases of the Kagome lattice Ising 
antiferromagnets with different kinds of spin dynamics. 
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any case, both the valence-bond crystal and spin liquid 
phases found here are gapped (gap energy scale may be 
very small), and the same is expected to be true for the 
isotropic Heisenberg model, regardless of what phase it 
actually lives in (unless it sits at the critical point). 

The Kagome lattice is a representative frustrated mag- 
net suitable for learning more general lessons on the two- 
dimensional systems. One important question, driven 
by efforts to discover unconventional Mott insulators, is 
under what circumstances can disordered and spin liq- 
uid phases be found in the frustrated spin models. One 
mechanism that clearly emerges is adding sufficiently 
strong further-neighbor and multiple-spin exchange pro- 
cesses. This has been already indicated in various other 
cases However, the corner-sharing lattices have 
been in focus due to a belief that even with the shortest- 
range dynamical processes one can still obtain the spin 
liquid physics. At least for the transverse field Ising 
model a disordered ground state is found in the Kagome 
system, making a sharp contrast to other usually stud- 
ied systems^ (most dimer models, triangular and fully 
frustrated square lattice Ising models). This disordered 
phase is conventional, since the system behaves almost as 
completely decoupled. Apparently, the completely local 
transverse field dynamics is unable to bring up correla- 
tions in a poorly connected lattice. However, as soon 
as the transverse field is replaced by the next least cor- 
related kind of dynamics (XXZ), a valence-bond ordered 
phase seems to emerge instead of a spin liquid. Since this 
happens in one of the most prominent systems for the 
short-range spin liquid, it is reasonable to speculate that 
the spin liquid in similar weakly perturbed quantum an- 
tiferromagnets quite generally requires further-neighbor, 
multiple-spin and ring exchange dynamics. Having a less 
connected lattice makes it easier for a trivial disordered 
phase to appear as a result of short-range dynamics, but 
perhaps not so much easier for the spin liquid phase. 

Also, the arguments from the end of the section IIV Fl 
indicate that the conservation of total Ising spin is a suf- 
ficient condition for non-trivial topology of disordered 
phases that have no net magnetic moment on any site. 
The SU(2) symmetric models, such as the Heisenberg 
model, are included as a special case, in agreement 
with the extension of the Lieb-Schultz-Mattis theorem 
to higher dimensions Clearly, the spin liquid can exist 
beyond this condition: as a stable phase, it can resist 
sufficiently weak spin non-conserving perturbations. 

On the technical side, this paper demonstrated an al- 
ternative U(l) gauge theory to that of the Ref.'^". Even 
though the present theory is more complicated, it pro- 
vides different insight, and avoids some difficulties that 
otherwise might have been encountered for the XXZ 
problem. One of its advantages is ability to give in- 
formation on the character of the ground (and excited) 
state wavefunctions, and a visual template for the kind 
of valence-bond orders that are possible in the Kagome 
XXZ models. We demonstrated a powerful analytical 
approach that successfully handles macroscopic degener- 



acy in the frustrated systems and extends the mean-field 
theories of the unfrustrated systems. 
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APPENDIX A: PROPERTIES OF THE TFIM 
LATTICE THEORY 

Here we derive two important properties of the lattice 
field theory H27fl with regard to its saddle-points. The 
saddle-point vectors ^ are given by H31|l , where the bond 
variables f^^^ describe the dimer coverings of the Kagome 
lattice with one dimer on every triangle, and an arbitrary 
even number of dimers on every hexagon. The value 
= 1 represents a dimer, while ^(y) = represents a 
vacancy. First, let us calculate the normalization of the 
saddle-point vectors ^ (the ^^^-^ variables have no time 
dependence, and we will drop the summation over time): 



(Al) 

31 J2 




jei 

const. + 4 



^42 + (.isCii 



In the last line we have used the notation from the 
Fig. [521 Switching to the bond variables, we have: 



oc const. - 4^(^^(ji)^(,2) + C(z3>?(j4>) (A2) 

i 

A \(«) / Hi) 



— const. 

— const. 



We have used the facts that the sum of ^^.y^ on every 
triangle is 1, since every triangle holds one dimer [N is 
the number of Kagome sites), and that the total number 
of dimers on the lattice is fixed (equal to the number of 
triangles). We see that all saddle-point vectors ^ have 
the same normalization. 

Now let us calculate how the coupling matrix C from 
the action acts on the saddle-point vectors ^. The 
quadratic parts of the expression H3U|) reveal how the 
matrix C acts on the height field vectors x whose com- 
ponents are Xi ^^nd A^. Substituting there X^jei ?u ^'^^ 
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FIG. 23: Local neighborhood of a site i; ei^ij) is —1 on the 
bonds of the shaded triangles, and +1 on the other bonds. 
Kagome bond orientations are also shown. 



every Xi, and f^'' every A,; reveals the action 

of C on the saddle-point vectors ^: 



(A3) 



where LiJ[ij) is neither a vector, nor a bond scalar: 



(A4) 



Since we need to find the sums of around a particular 
site, let us take a closer look at the neighborhood of a 



site i, and refer to the Fig. 1231 



fcej \ kej / 

4ey - 2(5j- + ^12) - 2(5j- 2(0z + 61) 
-2(5j,3(0* + C34) - ^SjAiiy, + ^43) 



(A5) 



Then: 



611^.^-2(62+61+64+^4 
6^6, ; (A6) 
^'Y^{i3)^ii - 2(e(ii)62 +e(i2)6i 
+e (i3> 64 + £(44)63) = ■ 



We see that the action of C on a vector whose compo- 
nents are Y.jei^v ^^"^ I]jei^{ii>^y simply reproduces 
those components, with an additional factor of 6. There- 
fore, all saddle-point vectors | are degenerate eigenvec- 
tors of the coupling matrix C, with an eigenvalue 6. 
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